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(N \ ABSTRACT 

Despite the increasing number of studies of barred galaxies at intermediate and high 
' redshifts, double-barred (S2B) systems have only been identified in the nearby (z ^ 

0.04) universe thus far. In this feasibility study we demonstrate that the detection and 
analysis of S2Bs is possible at intermediate redshifts (0.1 < z < 0.5) with the exquisite 
resolution of the Hubble Space Telescope Advanced Camera for Surveys (HST/ ACS). 
We identify barred galaxies in the HST/ ACS data of the Great Observatories Origins 
■ Deep Survey (GOODS) using a novel method. The radial profile of the Gini coefficient 

ly-j | - a model-independent structure parameter - is able to detect bars in early-type 

. galaxies that are large enough that they might host an inner bar of sufficient angular 

size. Using this method and subsequent examination with unsharp masks and ellipse 
fits we identified the two most distant S2Bs currently known (at redshifts z = 0.103 
and z — 0.148). We investigate the underlying stellar populations of these two galaxies 
through a detailed colour analysis, in order to demonstrate the analysis that could be 
,H , performed on a future sample of intermediate-redshift S2Bs. We also identify two 

S2Bs and five S2B candidates in the HST/ ACS data of the Cosmic Evolution Survey 
(COSMOS). Our detections of distant S2Bs show that deep surveys like GOODS and 
COSMOS have the potential to push the limit for S2B detection and analysis out by 
a factor of ten in redshift and look-back time (z w 0.5, At w 5Gyr) compared to the 
previously known S2Bs. This in turn would provide new insight into the formation of 



o 



Oh 



these objects 
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1 INTRODUCTION AND MOTIVATION 

The high redshift universe allows us to probe galaxy evolu- 
tion directly, by presenting us with examples of galaxies at 
a fraction of the age of the universe. The resolution needed 
to study these objects is possible either with ground-based 
adaptive optics or with the Hubble Space Telescope (HST). 
Using such observations, disk galaxies have been detected 
up to redshift z ~ 1, when the universe was not even half 
of its present age. How do the p roperties of these galaxies 
compare with those of local disks? [Lilly et all <ll998t) showed 
that the scale-lengths of the z ~ 1 disks were similar to those 
of local galaxies. Likewise, bulges at z ~ 0.5 are found to 
have had a small amount of luminosity evolution and no 



* E-mail: tlisker@astro.unibas.ch 
f Brooks Prize Fellow 



significant size evol ution assman et al.ll2002t . although 
ISimard et al.l 1119991) found a population of z ~ 0.5 galaxies 
which have higher bulge surface brightness than the local 
galaxies. 

Morphologies of high-redshift galaxies also have been 
studied and compared with those of local galaxies. The 
most prominent morphological fea ture of a large frac- 
tion of local galaxies is a bar ijS ellwoqd fc Wilk inson! 
Il993t lEskridee et"aH l2000l iKnauen et al l l2000tl . Since 
these can drive very strong secular evolutio n [e.g. 
iLvnden-Bell fc KalnaiJl972HDebattista et al1l200rl . know- 
ing how common they are at high z is important for con- 
straining disk evolution. Using a sample of 46 galaxie s ob- 
served in the Hubble Deep Field. lAbraham et al.Hl999t) con- 
cluded that the fraction of barred galaxies was ~ 24% at 
z ~ 0.2 — 0.7 but decreased sharply from there to z ~ 1.1. 
This result subsequently turned out to be an artifact of the 
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still limited resolution of the Wide Field Planetary Camera 2 
(WFPC2) , and several groups have now found a nearly con- 
stant barred fraction out to z ~ 1 using the higher resolution 
made possible by ACS (Elm egreen et all 120041 : Ijoeee et alJ 
12004^1 . Other disk structures which have been identi fied at 
high redshift include truncat ions at 0.6 < z < 1.0 iPerez 
| 2004 l :lTruiillo fc Pohlenl2005l) and warps dReshetnikov et al. 

Double-barred (S2B) galaxies, in which a small 'sec- 
ondary' bar is nested inside a larger 'primary' ba r, have 
inste ad been identified only out to z = 0.04 thus f ar fcrwinl 
2004). Indeed, the S2B catalog of lErwinl i2004) contains 
only two galaxies (out of 50) at a distance of more than 
100 Mpc and none further than 150 Mpc. S2Bs, which were 
first described by Ide Vaucouleursl lll975l) . have been pos- 
tulated to be a mechanism for driving gas past the inner 
Lindblad resonance of primary bars to feed the supermas- 
sive black holes (SMB Hs) that power active galactic nuclei 
JShlosman et alJll990D . Despite the possibility of their play- 
ing such an important role in SMBH nurture, understanding 
them has progressed significantly only in recent years. Sec- 
ondary bars have been sought in gas-poor early-type barred 
galaxies with high resolution HST imaging in order to dis- 
tinguish secondary bars from other nuclea r structures such 
as disks, rin gs, spirals and dust extinction ferwin fc Sparkd 
Il999l . 120021) . We now know that secondary bars are unam- 
biguously found in about 30% of early-type barred galax- 
ies. Moreover, we finally have direct observational proof 
for a kinematic decoupling between the pri mary and sec- 
ondary bars in one S2B galaxy, NGC 2950 (ICorsini et alJ 
2003). Such decoupling is necessary, but not sufficient, for 
secondary bars to funnel gas inwards. On the theoretical 
side, important progress has c ome f rom the development by 
iMacieiewski fc Sparkel Jl997l l200d) of the orbit formalism 
necessary for understanding their dynamics. Finally, simu- 
lators have s tarted to produce gas-f ree S2B systems in their 
simulations dRautiainen et ai1l2002l) . The secondary bars in 
these simulations rotated faster than the primary bars and 
survived for several gigayears, properties which favor SMBH 
fueling. 

Observing S2Bs at significantly earlier epochs would al- 
low us to compare them with local S2Bs and would thus 
provide new insight into the formation and evolution of 
these objects. Because of their small scales, their dynam- 
ical timescales may be quite short (~ 10 Myr). However, 
this does not necessarily imply a short lifetime: since as 
many as ~ 1/3 of early-type barred galaxies harbor sec - 
ondary bars (e.g. lErwin fc Sparkell2002t lLaine et al.ll200^ . 
it appears very unlikely that nuclear bars are destroyed and 
regenerated repeatedly in these nearly gas-free systems. In- 
stead, they probably have been in place a long time. It is 
thus worthwhile to investigate whether the relative proper- 
ties of the secondary and primary bars show any evidence 
for evolution. A comparison of the stellar populations of in- 
ner and outer bar would be one way to address this issue. If 
the stellar populations of secondary bars are systematically 
younger than those of the corresponding primary bars, this 
would support the hypothesis tha t nuclear bars form ou t of 
gas driven in by the primary bar JShlosman et al]ll989T) . If 
instead the secondary bars consisted of older stellar popula- 
tions than the primary bars, this would make this scenario 
implausible. 



Given that the difference between the optical proper- 
ties of two stellar population s decreases rapidly with age 
fe.g. iBruzual fc Charlotl f2003l. it would be desirable to ob- 
serve S2B galaxies at larger look-back times. Moreover, a 
large enough sample of S2Bs at higher redshifts would reveal 
a possible evolution of the ratios of secondary-to-primary 
sizes. If inner bars lose angular momentum to the outer bar, 
they might grow larger in radius, leading to a larger ratio as 
time goes on. 

While such questions motivate the search for and study 
of S2Bs beyond the nearby universe, it is not yet clear 
whether such a study would be at all possible with present- 
day instrumentation. Therefore, as a first step, it is neces- 
sary to perform a feasibility study, i.e. to show that a) S2Bs 
can be found at significantly larger look-back times than 
known to date, and b) the stellar populations of primary 
and secondary bars can be analyzed and distinguished. An 
estimate of the limiting redshift for detecting inner bars with 
th e HST Advanced C amera for Surveys (ACS) is provided 
bv lSheth et~ai] (120031) , who give a lower limit for a bar semi- 
major axis of 2.5 times the point-spread function (PSF) full- 
width at half- maximum (FWHM). Adopting 1 PSF FWHM 
= O'.'ll and the concordance cosmology (Ho = 70 km s _1 
Mpc - and Qa = 1 — J! m = 0.7, which we assume through- 
out), the resolution of the HST/ ACS would be sufficient to 
identify inner bars out to z ~ 0.25 provided t hat their size s 
are similar to those of local nuclear bars from lErwinl i2004l) . 
This limit is illustrated in the bottom panel of Fig. by 
the dotted vertical lines, and is significantly larger than the 
redshifts of known S2B systems, z ^ 0.04. 

Motivated by the foregoing estimate, we present a feasi- 
bility study for detecting and analyzing double-barred galax- 
ies at intermediate redshifts, including our identification of 
the two most distant S2Bs known so far. In Section|2]we de- 
scribe briefly the publicly available data used in this study. 
Section explains how the objects were found, including a 
novel technique for identifying bars. A structural analysis of 
primary and secondary bars is presented in Sectional fol- 
lowed by a colour and stellar population analysis in Section 
|S] Several (candidate-)S2Bs in the Cosmic Evolution Survey 
(COSMOS) are presented in Section|SJ and a final discussion 
is presented in Section |7] 



2 THE DATA 

The HST/ ACS images of the Great Observatories O rigins 
Deep Survey (GOODS) release vl.O jGiavalisco et al.l f2004l 
cover an area of 300 arcmin 2 divided into two fields: the Hub- 
ble Deep Field North (HDFN) and the Chandra Deep Field 
South (CDFS). They consist of four deep images through 
the F435W (B), F606W (V), F775W (i), and F850LP (2) 
passbands of ACS. The CDFS is divided into 18 sections and 
the HDFN into 17 sections with the images drizzled with a 
pixel scale of 30 mas, chosen to optimize the sampling of the 
PSF. The images are not perfectly background subtracted: 
there is a remaining background of the order of one third 
of the noise RMS, which can vary by up to a factor of two 
across the image. Since this can contaminate measurements 
of colour and morphological parameters, we subtracted an 
individual background flux value for each object and band. 
We defined this background as the median pixel value within 
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a box of 21" x 21" centered on the galaxy, after masking all 
sources. 



3 OBJECT SELECTION: GINI PROFILES 

Despite the common use of ellipse fits for bar identific ation , 
their limitations are well known Ce.g. lErwin fc Soarkelfc ofm 
and mainly originate from the fact that ellipse fits assume 
a definite (namely elliptical) isophotal shape of the object 
in question, which might not be the case in reality. Mo- 
tivated by the successful application of model-independent 
structural parameters such as asymmetry or concentration 
to ga laxy morphology (e.g. lAbraham et al.l|2003l : IConselicel 
2003), we sought a bar identification method that would 
rely solely on such parameters and would therefore provide 
a computationally quick and easy way to detect bars. 

Given the usual difficulties in reliably identifying nu- 
clear structures in late-type galaxies because of the signif- 
icant amount of patchy obscuration that could mimic bar- 
like central shapes, we restricted our study to early-type 
galaxies, i.e. E/SO galaxies as well as early-type spirals. 
We did this by selecting objects with central con centra- 
tion values C > 3.5 (cf. Fig. 11 of lLotz et alll2004) . where 
C = 51og[r(80%)/r(2 0%)l and r(x%) is the radius enclosing 
x% of the total flux feershadv et alJl200(t> . The total flux 
was measured in the i-band within a circular ape rture of 
r — 1.5x the Petrosian radius, Rp jPetrosianlll976h . where 
Rp is given by 



m(Rf 
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(1) 



and [i is the surface brightness. We set cutoffs at apparent 
magnitude m^AB ^ 22.0 and half-light radius ^ 10 pixels 
to ensure that the objects were bright and large enough for 
identification of an inner bar. These criteria yielded a work- 
ing sample of 231 objects in the GOODS dataset. 

To detect large bars with major axis exceeding 15 — 20% 
of the size of the galaxy, we computed model-independent 
radial profiles of structural parameters in circular apertures 
ranging from 1/30 Rp to 1.5 Rp on the i-band images. We 
focus on large bars for two reasons: first, the bar needs to 
stand out in the radial structure profile, hence requiring its 
size to be a significant fraction of the galaxy radius, and 
second, small bars are less promising candidates for hosting 
nuclear bars that are large enough to be detected. 

One structure paramet er - the Gini coefficient, G iGinil 
Il912l : lAbraham et al.ll2003F) - proved to be a useful diagnos- 
tic for larg e bars . We calculate G following the definition of 



tic tor large bar 
lLotz et alJ lEoO 



G = 
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(2) 



where N p i x is the number of pixels in the image and X is 
th e flux and Xi ^ X2 < X3 ^ ... < ^iv pix . As discussed 
bv lLotz et alJ (Hb04), the absolute values of X are required 
in order to preserve the correct structure at low flux levels, 
where noise can result in negative values of X after back- 
ground subtraction. G quantifies the distribution of flux val- 
ues among the pixels of an object's image. The two extreme 
values are G — if all pixels have the same flux, and G = 1 



if all the flux is concentrated in a single pixel. G can there- 
fore be thought of as measuring the contrast between the 
faintest and the brightest regions of a galaxy's image, no 
matter where in the galaxy they are located; thus it does not 
assume any shape for the object. In our definition G(r) is 
computed within each aperture as described above. Circular 
apertures are preferred over elliptical ones because the shape 
of a bar clearly stands out in a circular aperture. An ellip- 
tical aperture instead would fit the bar's shape to a higher 
degree, preventing the bar from standing out quite as much. 

What do the Gini profiles of ordinary galaxies look like? 
From the objects in our working sample we assembled a 
subsample of such galaxies, ranging from ellipticals to late- 
type spirals and irregulars based on a visual classification. 
Since we selected objects with high concentration values, 
late-type spirals appear in our working sample generally 
when they are seen nearly edge-on. Similarly, the objects 
classified 'irregular' still require a large amount of their light 
to be centrally concentrated to enter the sample. It is ap- 
parent from Fig. [5] that G(r) steadily increases with radius 
for E/S0 galaxies (possibly including edge-on early- type spi- 
rals), while for spirals it reaches a maximum and decreases 
thereafter. Late-type spirals can be distinguished because 
their maxima in G(r) are lower and are reached earlier than 
in early-type spirals 

In the Gini profile of several galaxies, a plateau is seen 
at intermediate radii (top left panel of Fig. We iden- 
tify such a feature by its mathematical properties: the sec- 
ond derivative changes from negative to positive, with the 
first derivative positive before and after, but becoming small, 
or possibly negative, in between. We approximate the first 
and second derivative by considering small (Rp/30) inter- 
vals, which are then smoothed by averaging over three data 
points at each step. An object was selected if its profile shows 
a plateau that is at least as strong as the one presented in 
the upper right panel of Fig. |5] This also requires the bar to 
be 'strong', i.e. a significant fraction of the total light needs 
to make up the bar shape. Bars that are so weak that they 
would hardly be seen without the help of unsharp masking 
or ellipse fits probably would not have enough effect on the 
Gini coefficient to produce a clear plateau in the profile. A 
plateau is evident in 46 objects of our working sample. Vi- 
sual inspection shows that 17 of them indeed have a large 
bar, 23 galaxies are unbarred, while for 6 galaxies we could 
not decide whether a bar was present based on visual in- 
spection. 

How can the shape of G(r) for a barred galaxy be un- 
derstood? First, the initial increase of G(r) with radius is 
caused by the brightest pixels being all located in the central 
region; thus increasing the radius adds only fainter pixels. 
This leads to a continuously increasing spread in the pixel 
flux values, and thus each Gini profile begins with increas- 
ing G(r), independent of galaxy type. However, a bar is then 
able to maintain a relatively high surface brightness out to 
a larger radius than an unbarred object, therefore the frac- 
tion of faint pixels remains lower, which probably leads to 
the plateau in G(r). Thus the Gini profile provides a simple 
way of identifying large, strong bars, and it can be computed 
easily and automatically for all sources detected in a deep 
survey such as GOODS. 

Although we did successfully apply this bar detection 
method to the GOODS dataset, we point out that it is not 
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yet optimized: the fraction of unbarred objects among the 
selected galaxies is ^ 50%, which would seem to not be 
a significant improvement over ellipse fits. Moreover, when 
visually inspecting all 231 galaxies of our working sample, 
11 objects possess a bar but were missed by our method, 
partly because their bar is relatively small, but also because 
of spirals and lopsidedness. A quick visual study revealed 
that the method could be improved further, for example 
by including a second parameter to exclude those systems 
which are asymmetric. Moreover, our above interpretation of 
the cause of the plateau in the profile might be too simplistic: 
it could well be possible that not (only) the bar itself is 
causing the plateau, but rather the configuration of a barred 
galaxy as a whole, i.e. the distribution of pixel flux values 
around the bar. A full study of these issues is beyond the 
scope of this paper and will be addressed elsewhere (Pasquali 
et al., in prep.). Since completeness is not the goal for our 
feasibility study, we are satisfied with the easy derivation of 
the Gini profiles and defer further quantitative evaluation of 
the method's potential. 

Can the method be applied recursively to identify nu- 
clear bars inside barred galaxies? This did not prove fea- 
sible due to the small apparent size of nuclear bars and, 
more importantly, the significant amount of smoothly dis- 
tributed underlying galaxy light which decreases the con- 
trast of nuclear bars (for the same reason, unsharp masks 
are such a useful tool for studying nuclear bars). Thus 
we had to inspect the selected barred objects visually and 
through unsharp masks in order to identify nuclear bars 
(Sect. |1J. Two objects were clearly identified as double- 
barred: GOODS J033230. 93-273923. 7, hereafter S2B 1, and 
GOODS J033233.46-274312.8, hereafter S2B 2; both lie in 
the CDFS part of the GOODS survey. 



4 STRUCTURAL PROPERTIES 

Here we present a structural analysis of S2B 1 and S2B 2. 
For th i s purp ose we use unsharp masks and ellipse fits as in 
lErwinl (1200411 . 



4.1 S2B 1 - a double-barred galaxy at z = 0.148 

The left panels of Fig.|3]show, from top to bottom, an i-band 
image of S2B 1, an unsharp mask, and the ellipse fits. The 
unsharp mask was created by dividing the original image 
with one convolved with a Gaussian of a — 40 pixels; in the 
center we insert another unsharp mask created through a 
convolution with a much smaller Gaussian (a = 4 pixels) to 
highlight the inner structure. Apart from the huge primary 
bar, a very prominent secondary bar is present along with 
a nuclear spiral. Its spectrosco pic redshift, as meas ured by 
the VIMOS VLT Deep Survey SLe Fevre et alJl2004l) . is z = 
0.148. S2B 1 is therefore the most distant S2B galaxy known, 
at a luminosity distance of 703 Mpc . The black scale bar 
in the image corresponds to a physical length of 5 kpc; the 
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image scale is 2.585 kpc/". We visually classify S2B 1 as 
SBO/a. 

In order to quantify the structure of S2B 1, we per- 
formed an ellipse fit using the 1RAF task ELLIPSE. The 
result is shown in Fig. |3] both the primary and secondary 
bar can be clearly identified as significant local peaks in the 
ellipticity. We did not attempt a deprojection of S2B 1, since 
its outer structure is difficult to determine: the primary bar 
dominates the galaxy's appearance, and it cannot be ascer- 
tained whether it is circular or still part of an oval distortion. 
Since the primary bar lies along the major axis of the outer 
isophotes, as indicated by the nearly constant position angle 
in the ellipse fits, we do not expect measurements of its size 
and shape to be significantly affected by a, probably modest, 
inclination. 

Following lErwinl J2004I) . we measured the primary bar's 
semi-major axis at its maximum ellipticity a e> i = 7.0 kpc, 
and the upper limit on the bar semi-major axis as Lbar.i = 
8.3 kpc, where Lt ar is the smaller of a m i n and aio, with a m i n 
being the semi-major axis at the first ellipticity minimum 
outside the bar, and aio being the semi-major axis where 
the position angle has changed by 10° from the bar posi- 
tion angle. The equivalent values for the secondary bar are 
a e ,2 = 0.8 kpc and Lb ar ,2 =1.1 kpc. The physical size of the 
primary b ar in S2B 1 i s larger than 47 of the 50 S2Bs pre- 
sented bv lErwinl ll2004l) . while the secondary is larger than 
76% of Erwin's sample. 

Visual inspection of the immediate environment of 
S2B 1 reveals a jet-like object close to it and pointing to- 
wards its center (Fig. . This object is seen in all bands. 
If the secondary bar in S2B 1 is feeding gas to an SMBH, 
then a jet originating from the SMBH may indeed be ex- 
pected. However, it s photometric red shift provided by the 
COMBO-17 survey JWolf et alJl2004h is z = 0.681, suggest- 
ing the object is a background galaxy. On the other hand, 
one would not expect a photometric redshift (which relies 
on template spectral energy distributions) to yield useful 
results if the object actually was a jet. The X-ray flux of 
the galaxy would allow us to test for an acti ve black hole 
in th e center. The Extended CDFS survey dVirani et alJ 
2006) does not detect S2B 1, but it provides an upper limit 
on the X-ray flux. The 3 sigma upper limit on the count 
rate is 6.35 x 10 counts/s (Virani, priv. comm.), which 
translates into an upper limit on the flux in the soft bands 
(0.5-2.0 keV) of 2.7xl0~ 16 erg cm -2 s -1 , using the Chandra 
X-ray Center Proposal Planning Toolkit PIMMS v3.7 2 with 
a power law spectral model and a photon index of 1.8. This 
corresponds to an upper limit on the soft X-ray luminosity 
of Lx = 1.6 x 10 40 ergs" 1 . B y applying equations 1 and 2 of 
iHornschemeier et alJ i2003T> we derive an upper limit on the 
X-ray-to-optical flux ratio of S2B 1 of log(/ x // R ) < -2.9. 
This value lies two orders of magnitude below the typica l 
values for AGNs (cf. Fig. 3 of IHornschemeier et al.l [2003): 
how ever, a heavily obscu red AGN cannot be excluded. 

lAfonso et alJ i2006T) presented Australia Telescope 
Compact Array radio (1.4 GHz) observations of the CDFS. 
They found a total of 64 radio sources within the HST/ ACS 
region of the GOODS/CDFS with flux densities between 
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63 /iJy and 20 mjy. However, S2B 1 was not detected in this 
survey. 

4.2 S2B 2 - a double-barred galaxy at z = 0.103 

The right panels of Fig. [3] show an i-band image, unsharp 
mask, and the ellipse fits for S2B 2, which we classify as 
SBO/a. The unsharp masks were produced with Gaussians 
of a — 1.5 pixels and a — 30 pixels. The small unsharp mask 
reveals an inner bar and spiral in the center of S2B 2. This is 
confirmed by a local peak in ellipticity, although this peak is 
relatively weak due to the small apparent size of the bar. The 
position angle nicely tracks the inner bar and then changes 
by about 90° outside of it, since the inner and outer bars in 
this galaxy are nearly perpendicula r. S2B 2's spectrosco pic 
redshift, as measured by 2dFGRS JColless et al.ll200ll) . is 
2 = 0.103, putting it at a luminosity distance of 475 Mpc. 
The black scale bar in the image corresponds to a physical 
length of 5kpc; the image scale is 1.893 kpc/". 

We deprojected the galaxy by assuming that the ellip- 
ticity minimum outside of the primary bar is the appar- 
ent ellipticity of the galaxy as a whole; we then obtained 
the deprojected size of the inner bar as a e .2 = 0.4 kpc and 
Ltar,2 = 0.5 kpc. The corresponding values for the primary 
bar are a £j i = 4.2 kpc and Lt ar .i ~ 6.4 kpc, which did not 
require deprojection since the primary bar is almost along 
the galaxy's apparent major axis. 

The primary bar of S2B 2 is slightly s maller than the 
median size of the 50 S2Bs in lErwirl J20041 . while the sec- 
ondary bar is smaller than that of 70% of the same sample. 
With a e , 2 = 0'.'17 = 1.5 PSF FWHM the apparent size of 
the secondary bar is obviously just large enough for it to 
be detected by HST/ ACS (1 PSF FWHM = O'.'ll in the i- 
band and similar in the other bands). We therefore conclude 
that z ~ 0.1 would be the maximum redshift for HST/ ACS 
detection of the inner bar for the smallest ~ 30% of the 
currently known S2Bs. 

S2B 2 has a faint (ps eudo) ring around the m ain galaxy, 
resembling the Ri class of lButa fc Crocker! il99ll) . Rings are 
quite common in barred galaxies, and simulations show that 
th ey can be produced at b ar resonances (see, e.g. the review 
bv lButa fc Combedll996h . 

S2B 2 is not d etected in the radio observations of 
lAfonso et all (§006), but is detected in the CDFS sur- 
vey o f the Chandra X-Ray Observatory ^Alexander et all 
120031: ICiacconi et"afl |2002>) . Its X-ray flux is 1.8 x 
10" 16 ergcm" 2 s _1 in the soft bands (0.5 - 2.0 keV), which 
translates into a soft X-ray luminosity of Lx = 4.9 x 
10 39 ergs -1 . The upper limit on the flux i n the hard bands 
(2.0 - 8.0keV) is 9.3 x 10" 16 erg cm" 2 s" 1 iAlexander et alJ 
2003). With an X-ray-to-optical flux ratio of log(/x//R.) = 
—3.4 S2B 2 belongs to the optically bright, X-ra y faint 
(OBXF) class of sources IIHornschemeier et al]l2003l) . These 
objects are thought to be quiescent, normal galaxies, with 
their high-energy emission originating from stellar processes 
or from low-level accretion onto SMBHs, but not from lumi- 
nous active galactic nuclei unless they are heavily obscured. 
Since the optical and X-ray positions only differ by 0'.'4 - 
which equals t he median positional u ncertainty of the faint 
X-ray sources \ Alexander et al]|2003ft - it is likely to orig- 
inate from the central region. While we cannot gain more 
insight into the processes generating the X-ray emission of 



S2B 2, more detections of intermediate-redshift S2Bs at X- 
ray wavelengths would test the hypothesis that gas in S2Bs 
is fed to SMBHs. 



5 STELLAR POPULATION ANALYSIS 

Comparisons of the available photometry from the ACS im- 
ages with population synthesis models allow us to estimate 
the age and metallicity distribution of the stellar compo- 
nent in the two galaxies. A pixel by pixel approach has the 
advantage of spatially resolving regions which may have un- 
dergone different histories of star formation. However, the 
uneven distribution of the signal-to-noise ratio (SNR) across 
the galaxy images results in complicated maps of uncer- 
tainties which are hard to interpret. In order to improve 
the spatially resolved analysis, we decided to perform an 
adaptive bi nning following the Optima l Voronoi Tessellation 
method of ICappellari fc Copinl (|2003() , which has already 
been applied successfull y to a study of early-type galaxies 
in GOODS/CDFS (iFerreras et al]l2005ft . A Voronoi tessel- 
lation divides a region starting from a distribution of points 
called generators. The area (or volume) is then binned in 
those regions closest to a given generat or. Motivated by the 
analys is of integral field spectroscopy, ICappellari fc Copinl 
( 2003) presented an algorithm optimised for the binning of a 
2D map into bins which yield a small scatter of the SNR with 
respect to some target value and which are kept reasonably 
round in order to preserve, as much as possible, the spatial 
resolution. When we apply this method to colour maps, we 
ensure that the subsequent stellar population analysis is as 
homogeneous as possible. 

We generated colour maps by registering the images and 
convolving them with the PSF corresponding to the other 
band used for a given colour - e.g. for a V — i colour map 
we convolve the i-band image with the l/-band PSF and 
vice versa. We measured the PSFs from a set of stars in the 
GOODS/CDF S field. A recen t analysis performed by the 
GEMS group jRix et alJl2004l) found that a universal PSF 
obtained from a set of bright but unsaturated stars is fully 
sufficient for all GOODS sections and for all positions inside 
each section. For the Optimal Voronoi Tessellation on the 
registered and convolved images, we chose a target SNR of 
20 in the ratio of the V to i fluxes, which amounts to an 
uncertainty of ±0.05 mag for the V — i colour in each bin. 
The same bins were used in all other colours (B — V and 
i — z); the average errors are 0.1 mag in B ~ V and 0.05 
mag in i — z. Note that S2B 1 does not have a B-band image 
because it falls close to the edge of the field. Figures and |S| 
show the colour maps for S2B 1 and S2B 2, respectively. The 
distribution of colours is somewhat patchier than expected 
from the above errors only (i.e. the RMS of the colour values 
exceeds the errors), particularly in S2B 1, which is partly due 
to dust and partly to population differences. 

Table [T] lists some integrated properties of our objects. 
The l/-band absolute luminosity was computed from the 
observed passband which gives the smallest K correction, 
i.e. F606W for both S2B 1 and S2B2, with a correction 
Kab ~ 0.04 mag. The spectral energy distribution for these 
corrections was obtained from a simple stellar population 
with solar metallicity, whose age best matches the observed 
colours. 
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5.1 S2B1 

A young stellar component in the inner spiral structure is 
readily visible in the V — i colour map in Fig. at around 
0'.'75 from the center. Ongoing star formation is also de- 
duced from a low-resolution spectrum of S2 B 1 from the VI- 
MOS VLT Deep Survey of the CDFS region l)Le Fevre et alJ 
l2004h . The spectrum spans the 5600-9000 A region and fea- 
tures Ha and [OIII] but no [Nil] emission, a signature 
of star-forming HII regions. Although these might coincide 
with the above region of young stars, no spatial resolution is 
possible due to the slit size (1" x 10") and to ground-based 
resolution. 

The age-metallicity degeneracy prevents us from getting 
an accurate estimate of the age and metallicity distribution. 
Nevertheless, we can test for gross differences between the 
stellar populations in the primary and secondary bars via 
a colour-colour diagram. The top panel of Fig. |7| shows the 
V — i versus i — z colour distribution of S2B 1. Each "pixel" 
corresponds to a Voronoi-tessellated bin with the targeted 
SNR of 20 in the ratio of V to i, which aims at a 0.05 mag 
photometric error in V — i. The actual tessellation produced 
an average error of 0.047 mag with a 0.033 mag RMS scatter. 
Regions in the inner bar are represented by filled dots. We 
selected those regions of the galaxy with an i-band surface 
brightness between 17.5 and 18.5 mag/arcsec 2 , so that we 
exclude the very inner region which appears to be affected 
by dust. The filled triangles give the colours in the outer 
bar, located in the region with surface brightness between 
21 and 22 mag/arcsec 2 . We select pixels in a 20° wedge 
aligned with the position angle of the outer bar. We also 
distinguish pixels in two radial bins: the light grey symbols 
correspond to distances r < l'.'S from the center while the 
dark grey symbols are pixels at r > l'.'8. 

For comparison, we show the colo urs of simple stellar 
popul ations from the latest models of iBruzual fc Charlod 
(2003). Each line represents an age sequence between 2 and 
10 Gyr, and three metallicities are considered, as labelled 
in the figure. The dashed lines are dustless populations, 
whereas the solid lines assume E(B — V) = 0.2 and a sim- 
ple dust reddening law dCharlot fc Falllfe oOO). The inner bar 
shows the bluest colours in i — z, which can be interpreted 
as indicating that its metallicity is lower than of the outer 
bar, and in addition either the stellar ages or the amount 
of dust is higher. However a unique interpretation is not 
possible with these data. The outer bar regions possess a 
larger colour spread than the inner bar, which presumably 
is caused mainly by dust (given the patchy distribution of 
colours discussed above), but may partly be due to radial 
variations in the stellar populations. 



5.2 S2B2 

This galaxy has a conspicuous ring of young stars at a pro- 
jected distance ~ 3.8 kpc (~ 2") from the center (Fig.lSJ. In- 
terior to this ring, the red colours are indicative of either old 
populations or dust extinction. The bottom panel of Fig. Q 
shows the colour-colour diagram of S2B 2. The inner region - 
shown as dots - corresponds to an i-band surface brightness 
between 17 and 18 mag/arcsec 2 . The filled triangles are pix- 
els in the outer bar - selected in a 20° wedge centered at the 
position angle of the bar with surface brightness between 20 



and 21 mag/arcsec 2 . Two radial bins are chosen in the same 
way as with S2B 1, and similar population synthesis model 
predictions are shown. Note the remarkable difference be- 
tween the outer bar and the red colours of the inner bar, 
which is significantly affected by dust. The figure hints at 
younger ages in the outer parts of the primary bar, which is 
clearly due to the ring of young stars mentioned above. 



6 S2BS IN COSMOS 

In the course of our study, the reduced HST/ ACS i-band 
imaging data (vl.2) of the first square degree of the Cosmic 
Evolution Survey (COSMOS 3 ) became publicly available. 
Since these data cover an area about twelve times larger than 
GOODS, they can be expected to host several S2B galax- 
ies at intermediate redshifts. The COSMOS image depth is 
only slightly lower than that of GOODS, therefore a search 
for S2Bs with our Gini profile method and subsequent vi- 
sual inspection appeared promising. We applied the same 
selection criteria as for our GOODS dataset, i.e. C ^ 3.5, 
7Yii ] ab ^ 22.0, and half-light radius ^ 10 pixels. Indeed, 
we were able to identify two S2Bs and five S2B candidates 
which we list in Table H Their images and unsharp masks 
are shown in Fig. |H| 

Based on their unsharp masks and ellipse fits (not 
shown), the objects COSMOS J095922.76+024245.3 and 
COSMOS J095936.98+015107.6 clearly host inner bars. 
Four of the other five galaxies probably contain nuclear disks 
since the inner position angles match those of their outer 
disks. The remaining object is an S2B candidate. In this 
analysis, the larger pixel scale of the COSMOS images as 
compared to GOODS has a noticeable effect - a drizzling 
process that would bring the COSMOS pixel scale from the 
current 50 mas down to 30 mas as in GOODS would re- 
sult in a better discrimination of inner bars. Unfortunately, 
no further data (such as colours and redshifts) are publicly 
available for these galaxies; therefore their distances, bar 
sizes, and stellar population properties remain unknown. 



7 SUMMARY AND DISCUSSION 

To date, all observational studies of S2Bs have been under- 
taken in the nearby universe. In this paper we presented 
the two most distant S2Bs currently known, at z = 0.103 
and z = 0.148, corresponding to a look-back time of 1.3 
and 1.9 Gyr, respectively. Although our S2B sample con- 
tains only two objects, we can perform a consistency check 
with the number of S2Bs that would be expected given the 
numbers for nearby objects. Our GOODS sample of early- 
type barred galaxies consists of 17 objects that were se- 
lected by our Gini profile method, plus 11 barred galaxies 
that were missed by the method (see S ect. El. Adopting an 
S2B fraction of 30% (e.g. lErwin fc SparkeN2002ft we would 
thus expect 8.4 S2Bs to be within our sample. With a given 
distribution of redshifts, the bottom panel of Fig. then 
leads to an estimate of the number of S2Bs whose inner bar 
would be large enough to be identified. Redshifts are avail- 
able for 12 SBs in GOODS/CDFS from the COMBO-17 

3 See http://www.astro.caltech.edu/cosmos 
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survey J Wolf et alJl2004) and the VIMOS VLT Deep Sur- 
vey (VVDS. iLe Fevre et all 12004) . with 6 objects beyond 
z — 0.5. When we adopt this distribution to be representa- 
tive for all 28 SB galaxies, we would expect to detect 2.4±1.5 
S2Bs, in excellent accordance with our findings. While this 
suggests no significant evolution of the S2B fraction, clearly 
a larger sample of S2Bs at intermediate redshifts is necessary 
to address the statistics properly. 

Any sample of intermediate-redshift S2Bs would of 
course be biased towards large secondary bars, and also to- 
wards large primary bars - the latter being due to both 
the selection of barred galaxy candidates and due to a cor- 
relation of the sizes of inner and outer bar. Despite this 
inevitable bias, it would be possible to compare the rela- 
tions between the sizes of primary and secondary ba r with 
those of the local sample: the data from lErwinl (|2004) show 
a clear correlation of primary and secondary bar size (left 
panel of Fig. Moreover, the secondary-to-primary ratio 
is not constant but it increases with increasing secondary 
bar size (right panel of Fig. |§J . If an intermediate-redshift 
sample showed an offset with respect to the local relations, 
this would hint at a decoupled evolution of sizes, possibly 
related to the loss of angular momentum. While S2B 2 (open 
square in Fig. ^ falls onto the above relations, the primary 
bar of S2B 1 (open triangle) is very large compared with its 
secondary. Still, S2B 1 falls within the scatter of the rela- 
tions. 

In addition to the structural parameters, we presented 
a colour analysis of the two S2Bs. Unfortunately, it is not 
possible from the available data alone to precisely trace back 
the star formation history of an individual galaxy (and of 
its bars), especially given the possible effect of dust and 
the well-known age-metallicity degeneracy. Nevertheless, we 
show that population differences between the outer bar and 
the inner bar region can be quantified, demonstrating the 
feasibility of stellar population analyses of S2Bs at interme- 
diate redshifts. It is not clear whether the evolution of the 
primary bar and the secondary's stellar content are linked in 
some way. For example, the secondary bar could be formed 
from stars born from gas swept by the primary bar. An 
intermediate-redshift S2B sample of significant size would 
be able to address this question. If the amount of dust were 
small, the colours of S2B 1 could be interpreted as the inner 
bar having lower metallicity and larger stellar ages (Sect. 
15. II . making the above scenario implausible for this galaxy. 
The colours of S2B 2 are clearly more strongly affected by 
dust and little can be said about the relative ages of the 
primary and secondary bars. 

Another important result for future searches of S2Bs 
stems from our detection of S2B 2: its inner bar has an ap- 
parent semi-major axis size of only a Ei 2 =1.5 PSF FWHM 
but can still be recognized. This is significantly smaller 
than the estimate of 2.5 PSF FWHM as a lower limit 
from lSheth et all (12003ft . probably also due to the increased 
resolution from image drizzling. We thus suggest 1.5 PSF 
FWHM as a more modest minimum size for detecting nu- 
clear bars. When we adopt this new limit, Fig. Q demon- 
strates that we increase the maximum redshift for detection 
of the largest inner bars (a e ,2 > 1 kpc) of local S2Bs by more 
than a factor of two, from a redshift of z fa 0.25 to z ~ 0.5. 
In particular, S2B 1 with a e , 2 = 0'.'3 = 2.7 PSF FWHM 
should still be detectable at z ~ 0.3 given the new detec- 



tion limit. To test this, we created a mock image of S2B 1 at 
z — 0.3. Taking into account the corresponding wavelength 
shift would require an image right between U-band and i- 
band to be 'redshifted'. We therefore performed the artificial 
redshifting for both V- and i-band images. The images were 
rescaled and convolved with a Moffat profile representing 
the PSF of the GOODS images. Gaussian noise correspond- 
ing to the noise of the original images was then added, and 
unsharp masks, ellipse fits and Gini profiles were produced. 
The flux values were rescaled correctly throughout the pro- 
cess, taking into account surface brightness dimming as well 
as a proper normalization of the convolved images. We did 
not deconvolve the original image prior to redshifting it, 
which means that our final PSF is slightly too large; hence 
our treatment is conservative. The secondary bar can still 
be clearly identified in the unsharp masks (Fig. 1101 , and the 
primary bar would also have been selected from its Gini pro- 
file (dotted line in top left panel of Fig. The ellipse fits 
(Fig. llOt show the inner bar as a flat peak in ellipticity, with 
a sharp drop outside of it. This inner part looks very similar 
to the ellipse fit of S2B 2, for the obvious reason that both 
inner bars in the redshifted S2B 1 and in S2B 2 are close 
to their detection limit. Since the inner bar of S2B 1 is not 
excepti onal compared to local S2Bs (24% of the inner bars 
listed in lErwirJ j2004h are even larger), those could therefore 
be detected out to still larger redshifts, z ~ 0.5 (Fig. 0. 

Our identification of two S2Bs and five more S2B- 
candidate objects in the COSMOS data proves that in 
present-day deep surveys, it is worthwhile to search for dis- 
tant S2Bs. These galaxies can be analyzed in a similar way 
as presented in this paper as soon as redshifts and colour 
data become available. We have shown that HST/ ACS data 
are capable of making considerable progress in understand- 
ing the structure and evolution of such objects. We therefore 
see the possibility of an S2B sample distributed over a large 
redshift range in the near future. Such a sample would serve 
as an important observational constraint on models of S2B 
formation. 

We described a first step in developing a model- 
independent method for detecting large and strong bars 
which we used here. While the method as presented here 
is not optimized, it appears that with further refinements 
it can become a useful method for identifying bars in large 
surveys. 
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Table 1. Photometric properties of S2B 1 and S2B 2. All magnitudes and colours are given 
in the AB system. Magnitude errors are estimated to be 0™1. Formal errors from the SNR arc 
much smaller. Similarly, formal errors on the colour values are negligible. Instead, we refer the 
reader to Figs. [S] and [S] concerning the colour variation within the galaxies. See text for further 
details. 

Object ID ms M v B - V V - % i~ z 

magAB magAB mag A B mag A B mag A B 

GOODS J033230. 93-273923. 7 S2B 1 17.9 -20.8 0.57 0.27 

GOODS J033233. 46-274312. 8 S2B 2 17.1 -20.6 1.31 0.57 0.28 



Table 2. S2Bs and S2B candidates in COSMOS. Magnitude errors are estimated to be 
DPI. 



Object 


m 


Classification 




magAB 




COSMOS J095922. 76+024245.3 


18.4 


S2B 


COSMOS J095936. 98+015107.6 


17.9 


S2B 


COSMOS J095931. 25+015916.5 


18.1 a 


possibly S2B or SB with nuclear disk 


COSMOS J095947.65+021549.0 


18.4 


possibly S2B or SB with nuclear disk 


COSMOS J100045. 42+021932.4 


19.4 


possibly S2B or SB with nuclear disk 


COSMOS J100048.23+015856.2 


20.0 


possibly S2B or SB with nuclear disk 


COSMOS J095952. 25+020543.1 


20.1 


possibly S2B 



a Partially contaminated by a bright star's bleeding trails. 
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Figure 1. Bar sizes. Comparison of deprojected bar sizes (a s ) of the 50 nearby S2Bs from lErwinl J2004) to our objects. Top: Primary 
bar sizes. Bottom: Secondary bar sizes. Thick vertical lines denote the values for our double-bar detections (S2B 1 and S2B2). Vertical 
lines with redshift labels give the minimum physical size of a secondary bar that could still be identified at the respective redshift in 
the GOODS images or images with similar characteristics (pixel scale 30 mas, 1 PSF FWHM = O'/ll). The dot ted vertical lines assume 
a minimum angular size of the secondary bar's semi-major axis of 0'.'28 or 2.5 PSF FWHM iSheth et alJl2003l) . while the solid vertical 
lines assume a minimum angular size of 0'.'17 or 1.5 PSF FWHM, which is the semi-major axis of the smallest secondary bar we detect 
(S2B2). 
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Figure 2. Gini profiles. The radial variation of the Gini coefficient (G) is shown for our two S2Bs (upper left panel), as well as for 
a sample of unbarred galaxies (middle and bottom panels) which are visually assigned to different classes as indicated in each panel. 
The dotted line in the upper left panel is for S2B i after artificially redshifting it to z = 0.3 (see Sect.|7J- Arrows indicate a plateau in 
the profile (see text for details). The upper right panel shows an early-type barred galaxy whose profile has a weak plateau just strong 
enough to be selected. The uppermost profile in the middle right panel belongs to an early-type spiral seen almost egde-on. 
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Figure 3. Images, unsharp masks & ellipse fits. For S2B 1 (left panels) and S2B 2 we present i-band images (top panels), unsharp 
masks (middle panels), and ellipse fits (lower panels). In the images, north is up and east is to the left; position angle is measured cast 
from north. The black scale bar corresponds to a physical length of 5 kpc. Unsharp masks have been computed by dividing the original 
i-band image by a Gaussian-convolved image. Two unsharp masks are superposed: one using a large Gaussian <r (<r = 40 pix for S2B 1 
and cr = 30 pix for S2B 2), and a second one using a small value (cr = 4 pix for S2B 1 and a = 1.5 pix for S2B 2) which is shown as inset 
in the image center and is enlarged in the lower right corner for S2B 2. Ellipse fits (using the IRAF task ELLIPSE) are shown in black 
for the i-band images, with errors as given by ELLIPSE. The other bands are only shown as lines: dash-dotted for B-band, dotted for 
V-band, and dashed for z-band. The dotted vertical line indicates the value of 1 PSF FWHM (0"11). 
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Figure 4. A jet? i-band image of S2B 1 and the jet-like object pointing back to its center. The black scale bar corresponds to a physical 
length of 5 kpc. 
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Figure 5. Colour maps of S2B 1. An adaptive binning based on the Optimal Voronoi Tessellation of ICappellari fc Coninl l200.j) has 
been applied. A target SNR = 20 was chosen in the flux ratio between the V and the i passbands. The contours give the i-band surface 
brightness from 18 mag/arcsec 2 (inner contour) in steps of 0.5 mag. 
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Figure 6. Colour maps of S2B 2. Similar to Figure 151 
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V-j 

Figure 7. Colour-Colour diagrams. Colours are given separately for three regions in each galaxy (see 6" X 6" inset). Pixels in the 
outer bar are selected from a 20° wedge aligned with the bar and chosen to avoid dust lanes. Filled triangles represent regions in the 
outer bar at radial distances r < l'.'S (light gray, pointing upward) or r > T.'8 (dark gray, pointing downward). The filled dots are the 
colours at the position of the inner bar, chosen between 17.5 and 18.5 mag/arcsec (S2B 1) and between 17 and 18 mag/arcsec 2 (S2B 2). 
We exclude the very central region which could be significantly reddened by dust. The dashed and solid lines track the colours of simple 
stellar populations with E(B — V) = and 0.2, respectively. Each line is an age sequence from 2 to 10 Gyr. See text for more details. 
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Figure 8. S2Bs and S2B candidates in COSMOS. Objects are shown in the same order as in Table l2j i.e. the top two objects are 
classified as S2B, the middle four objects as possible S2B or SB with nuclear disk, and the lowermost object as possible S2B. For each 
object, the i-image is shown on the left and the unsharp mask created with a Gaussian of <r = 1.5 pixels on the right. North is up and 
east is to the left. The black scale bar corresponds to 1". The masks only cover the central area of each object, as indicated by the box 
edges in each image. 
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Figure 9. Size relations of inner and outer bars. Comparison of deprojected sizes (a e ) of inner and outer bars for the 50 nearby 
S2Bs from Erwinl I2004L filled dots) and for S2B 1 (open triangle) and S2B 2 (open square). Left: Primary versus secondary bar size. 
Right: Ratio of secondary-to-primary size versus secondary bar size. 
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Figure 10. Simulated observation of S2B 1 at z = 0.3. Unsharp mask (<j = 1.5 pix) and ellipse fits of S2B 1 after artificially 
redshifting it to 2 = 0.3. Both V-band (upper mask image; dotted grey line) and i-band (lower mask image; solid line) are shown. The 
dotted vertical line indicates the value of 1 PSF FWHM (O'.'ll). The corresponding Gini profile is included in Figure |2] See text for 
details. 



